arXiv:hep-ph/9507284v2 15 Nov 1995 


UTPT-95-11 

CMU-HEP95-06 

DOE-ER/40682-96 

JHU-TIPAC-950007 

hep-ph/9507284 

Hadron Spectra for Semileptonic Heavy Quark Decay 


Adam F. Falk 

Department of Physics and Astronomy, The Johns Hopkins University 
3400 North Charles Street, Baltimore, Maryland 21218 U.S.A. 

f alkSplanck.pha.jhu.edu 

Michael Luke 

Department of Physics, University of Toronto 
60 St. George Street, Toronto, Ontario, Canada M5S 1A7 

lukeSmedb.physics.utoronto.ca 


Martin J. Savage 

Department of Physics, Carnegie Mellon University 
Pittsburgh, Pennsylvania 15213 U.S.A. 
savage@thepub.phys.emu.edu 
(July 1995) 

Abstract 


We calculate the leading perturbative and power corrections to the hadronic 
invariant mass and energy spectra in semileptonic heavy hadron decays. We 
apply our results to the B system. Moments of the invariant mass spec¬ 
trum, which vanish in the parton model, probe gluon bremsstrahlung and 
nonperturbative effects. Combining our results with recent data on B me¬ 
son branching ratios, we obtain a lower bound A > 410 MeV and an upper 
bound m)( ole < 4.89 GeV. The Brodsky-Lepage-Mackenzie scale setting proce¬ 
dure suggests that higher order perturbative corrections to the first moment 
of the hadronic invariant mass spectrum are small for bottom decay, and even 
tractable for charm decay. 
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I. INTRODUCTION 


Our understanding of inclusive decays of hadrons containing at least one heavy quark 
has improved greatly over the last few years. The energy released during semileptonic or 
radiative decay of heavy hadrons is much larger than the scale Aqcd of the strong inter¬ 
actions, and therefore an operator product expansion (OPE) exists for some observables in 
these decays, including rates and differential spectra [il],||. The leading power corrections to 
the rates and lepton differential spectra for semileptonic decays of heavy hadrons ||-||] have 
been studied extensively, as have the power corrections to radiative decays 0. 

A major result of this analysis is that, except in regions where the expansion becomes 
singular such as the endpoint of the electron spectrum in semileptonic b —> u decay, the 
corrections to the parton model are quite small, suppressed by 0(Aq C - D /ml). While this 
does mean that the parton model is quite successful, it makes it difficult to test quantitatively 
the corrections given by the OPE. In particular, if quark/hadron duality does not hold in 
this energy regime, one would expect to see corrections to the parton model which could not 
be accounted for by the leading perturbative and 1/rriQ corrections. Shifman has recently 
criticized the related OPE analysis of r decay on the basis that violations of duality in the 
Minkowski regime introduce large corrections which are not seen at any finite order in the 


OPE g. 

In this paper we suggest that hadronic variables, in particular moments of the invariant 
mass spectrum dT/ds# and the hadron energy spectrum dT/dE H , provide a useful testing 
ground for the OPE. This is similar to the analogous suggestion, and analysis, for semilep¬ 
tonic r decays ||.|I(|. However, unlike the case for r decays, at tree level the final hadronic 
state at the parton level consists of a single quark. Therefore at lowest order in the OPE the 
final hadronic state has fixed invariant mass sh = m 2 q , and positive moments of (s# — m//), 
which are calculable as a double expansion in a s (mb ) and 1 /mb, directly probe physics be¬ 
yond the parton model. Similarly, at leading order in the OPE the maximum hadron energy 
is (ml +m^)/2mi, (when the quark q recoils back-to-back with the leptons); the region above 
this endpoint is populated only by gluon bremsstrahlung and nonperturbative effects. 

In this paper we calculate the corrections to the parton model results for these ob¬ 
servables, up to 0(l/ml,a s /mb). As discussed in Ref. [|lTJ , although the leading power 
corrections to leptonic variables arise at 0(1/ml), for kinematic reasons the leading power 
corrections to moments of the invariant mass spectrum arise at 0(1/mb). The 0(1/ml) 


corrections to the differential hadronic energy spectrum were first examined in Ref. 
however, we disagree with the results presented in that work. We also use the results of 
Ref. 0, in which the one-loop corrections to the hadron energy spectrum were calculated. 


We combine our results with recent data on B meson branching ratios to obtain a lower 
bound on the nonperturbative parameter A, which is the leading contribution to the differ¬ 
ence between heavy quark and heavy meson masses. 

Finally, using the BLM prescription |T4| to estimate the size of the two-loop perturbative 
corrections to the moments of the invariant mass spectrum, we demonstrate that the first 
moment appears to have a well-behaved perturbative expansion not only for B decays, but 
also for D decays, when the results are expressed in terms of physical observables. This 
suggests that studying hadronic observables in semileptonic decays of charmed hadrons, 
which are dominated by only two or three resonances, may shed some insight into the 
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FIG. 1. The kinematics for B —> X q lv. 


applicability of quark/hadron duality at low energies. 

II. KINEMATICS 

We start by introducing the kinematic variables describing the final state hadrons. For 
definiteness, we will consider semileptonic B decay, although the analysis extends simply to 
the decays of charmed hadrons. 

The kinematics of the inclusive process B —> X q £v is shown in Fig. [Tj. The four- 
momentum of the B meson is P£ = and q M is the four-momentum of the lepton 

pair. We write the four-momentum of the b quark as P{f, and assign the heavy quark the 
same four-velocity v fl as the heavy meson. The total energy of the leptons in the B rest 
frame is v ■ q, and their invariant mass is q 2 . It is convenient to define dimensionless parton 
level quantities Eq and So, 


E 0 = v ■ (P b - q)/m b = 1 - v ■ q , ( 2 . 1 ) 

so = (Pb ~ qf/ml = 1-2 vq + q 2 , 

where q M /m b . At leading order in 1 /m b , E 0 and s 0 are simply the scaled energy and 

squared invariant mass of the final hadronic state. However, since they are scaled by the b 
quark mass, this identification does not hold at subleading order in l/m b .f] Instead, they are 
related to the physical hadronic energy and squared invariant mass, 

E h = v ■ (. P B ~ q) = m B - v ■ q , (2.2) 

Sh = (Pb ~ q ) 2 = rrt 2 B - 2 m B v ■ q + q 2 , 


through 


— Ai T 3 A 2 ( x Ai T 3A 2 \ x, 

E h = A---1- I ms — A -|---] E 0 + ... , 

Zms \ ZrriB ) 

s H = m 2 + A 2 + (m| - 2km B + A 2 + Ai + 3A 2 ) (s 0 - rh 2 ) 
+ (2km b — 2A 2 — Ai — 3X 2 )E 0 + ... , 


(2.3) 


1 This fact was neglected in Ref. [|li|. 
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where rh q = m q /m b , and the ellipses denote terms higher order in 1 /m b . The quantities A, 
Ai and A 2 arise in the relationship between the quark and meson masses [|l5l , |T6il , 


- Ai + 3A 2 

m B = m b + A---h • • • , 

2m b 

t Ai — A 2 

m B * = m b + A--- 

2m b 


(2.4) 


From the measured B-B* mass splitting, A 2 ~ 0.12 GeV 2 . We note that in contrast to 
the lepton spectra, there are 1 /ra& corrections both to the physical hadronic invariant mass 
spectrum and to the physical hadronic energy spectrum, although these corrections are 
absent for the §q and E 0 spectra 
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While the complete shape of the E 0 spectrum may be calculated (away from the parton 
model endpoint E 0 = |(l + m 2 )) with the standard OPE analysis, only suitably averaged 
features, such as moments, of the So spectrum may be computed reliably. The difference 
arises because each point of the Eq spectrum receives contributions from states of different 
invariant masses, making the process inclusive, whereas by definition each point of the s 0 
spectrum only receives contributions from states of a single invariant mass. This may be 
seen explicitly by carrying out the usual OPE analysis for inclusive decays in the variables 
so and Eo , instead of the usual leptonic variables v ■ q and q 2 . 

The inclusive B meson decay rate is given by 


r (B -> X,lv) ~ J dS„ d£„ 0. iyiunso, 7 ), 


(2.5) 


where L^ v is the spin summed lepton tensor v oc {q^q u — Q^q 2 )- Using the optical theorem, 
the nonperturbative hadronic tensor is related to the imaginary part of the forward 
scattering amplitude [3],[2|, 


W" = y (B| |X)(X| r h \B)(2^)‘5 1 (P b - P x - q) 

X 


= —2 Im (B | i 
= 2 Irn T^ w 


d xe- iqx T 


jf(x)r h (o) I B) 


( 2 . 6 ) 


where — 75 )b. The time-ordered product T^ u may be written via an operator 

product expansion as a power series in a s (mb ) and 1 /m b . 

In the v ■ q plane, for fixed q 2 , the correlator T^ u has the analytic structure shown 
in Fig. |2]a, as discussed in Ref. j[||. There are cuts along the real axis, a physical one 
(corresponding to B decays) for v ■ q < \ [1 + q 2 — m 2 ], and an unphysical one (corresponding 
to scattering processes) for v ■ q > |[(2 + m q ) 2 — q 2 — 1]. The one-particle pole lies at 
the right hand end of the physical cut. After an integral over the charged lepton energy, 
the decay rate is computed by performing an integration over the top of the physical cut, 
for y/q 2 < v ■ q < |(1 + q 2 — m 2 ), followed by an integration over 0 < q 2 < (1 — rh q ) 2 . 
Note that in the limit rh q —> 0 and q 2 —> 1, the physical and unphysical cuts pinch the 
region of integration. In this corner of the parameter space, the operator product expansion 
breaks down. Attempts to resum the OPE to all orders in this region have thus far proven 
inconclusive P| JTTjJT8l . 
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(a) 


(b) 


physical cut 

unphysical cut 

one-particle pole 
I—►-1 

integration contour 


FIG. 2. The analytic structure of , (a) in the v ■ q plane, with q 2 fixed, and (b) in the §q 
plane, with Eq fixed. Both the physical and unphysical cuts are shown, as well as the position of 
the one-particle pole. 


Mapping from the v-q plane to the So plane at fixed Eq, one finds two cuts on the positive 
So axis, as shown in Fig. |2b. The physical cut, which terminates in the one-particle pole, 
extends over s 0 > m 2 q . The unphysical cut lies away from the pole, at so > rh 2 + 4 rh q + 4 Eq. 
The region of integration in so is given by max (2 Eq — 1 ,rh q ) < §o < Eq, to be followed by 
an integration in Eq, over rh q < Eq < 1. 

In these variables, the region of integration only touches the unphysical cut in the limit 
m q —> 0 and Eq 


0, which from Eq. (|2.1|) is equivalent to the condition for the cuts to pinch 
in the v ■ q plane. In this singular region, as before, the operator product expansion breaks 
down. We also note that the integration region covers the one-particle pole at So = only 
if Eq < -1(1 + m 2 ). Indeed, this corresponds to the maximum energy the final quark can 
take away in the decay process. The cut for s 0 > r ' n ‘q is populated only by multiparticle final 
states generated by the radiation of gluons. In perturbation theory, then, the differential 
spectrum dT/d E 0 for Eq > |(1 + rh 2 q ) is of order a s . 

As is the case for r decays, the contour of integration in Eq. Q2.5p may be deformed away 
from the physical region, except at the point the contour crosses the physical cut. However, 
we note that in contrast to r decays, the integrand in Eq. (|2.5|) does not have a double zero 
where the deformed contour approaches the physical region. It is possible, therefore, that 
deviations from quark/hadron duality in the Minkowski regime may be more pronounced in 
semileptonic heavy hadron decay than in r decay. 


III. SPECTRAL MOMENTS 

In this section we compute the spectral moments at the parton level. We will treat 
both the leading power corrections, proportional to Ai and A 2 , and the leading perturbative 
contributions, proportional to a s ( r m b ). We take the two types of corrections in turn. 
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A. Power Corrections 


For the computation of the power corrections, it is convenient to decompose the time 
ordered product into the form factors 

T^(s 0 , E 0 ) = -g^T.iso, E 0 ) + uW 2 (s 0 , £<,) + •■■, (3-1) 


where the omitted form factors vanish for massless leptons 
and T 2 the differential spectrum is given by 




dr 


To dsod-Eo 
where 


= lm\lEl - s 0 

7T 


\/ Eq — So 3(1 — 2Eq + so)7i(so, Eq) + {Eq — sq)T2(sq, Eq) , (3.2) 


r n = 


GW, \V bq \- 
192vr 3 


(3.3) 


is proportional to the total decay rate. 

The leading 1 /m^ corrections to the hadronic quantities T\ and T 2 were calculated in 
Refs. |3|,[|]. In terms of s 0 and Eq, they are given by 


Ti(sq, Eq) — 


s 0 - m 2 + ie 


Eq 

T 


Ai 

12m? 


Ad 


4 m? 


+ 


+ 


(s 0 - m 2 q + ie 


Al (5E 0 2 - 3E 0 - 2s 0 ) + A2 


T2(Sq, Eq) — 


(•So 

; 


m q + ie 


6 ml 

2Ai 

3m? 


2 ml 


{^Eq + Eq — 2sq) 


Eq(Sq — Eq) 


s 0 - rh 2 q + ie 


5Ai 5A 2 
6m 2 


+ 


+ 


(s 0 -m 2 q + ie ) 


(s Q - m 2 + ie 


'b 

7A, 

3m 2 

4Ai 

3m? 


2?n 2 


(^o — 1) H-^(5^0 — 3) 

mi 


(So - ^o 2 ) 


(3.4) 


(3.5) 


Integrating this expression with respect to Eq, we hnd the leading power correction to the 
invariant mass spectrum. Of course, since there is only a single quark in the final state, 
this expression is a singular function with support only at s 0 — m 2 . Only its moments, 
which we present below, are meaningful. The corrections to the hadronic energy spectrum, 
obtained by integrating first with respect to s o, are more interesting, and are presented in 
Appendix A. 

Because the expansions of T\ and T 2 in terms of l/m& contain pole factors l/(s 0 — m q ) n , 
it is simplest to compute the moments of (s 0 — m q ) rather than those of Sq. The requisite 
calculations are straightforward but tedious, and we present only the final results. It is 
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FIG. 3. The Feynman diagrams which contribute to the moments at order a s . There are also 
wave function corrections, which we do not show. 


convenient to scale the various contributions to T 0 rather than to the full width T ; the 
quantities which we will present below are then of the form 


J\/[ (n,m) 


1 

ro 



dr 

ds 0 dE 0 


ds 0 dE 0 , 


(3.6) 


for integers n and m. They are related to the parton level moments by a scaling to the 
corrected decay rate, 


(K(§o 


777 , 


;)”) = 


n,m) 


(3.7) 


B. Perturbative Corrections 


The perturbative corrections to T\ and T 2 are most conveniently calculated directly 
from the graphs in Fig. |3j. The radiative contributions to ((sq — m^) n ) come only from 
bremsstrahlung graphs and are straightforward to compute for arbitrary rh q . We find 


dT 
T 0 dso 


0‘s 1 [ (*,- 1 ) 

7T s 0 - m 2 q [ 27s§ 



6 Til q + s 0 (18?hg + 81r7 ig + 48?hg) 


(3.8) 


+ Sq(93 - 316mJ + 243 m\ + 102m®) + sjj(-41 - 478mJ + 9mj) 
+ Sq(—95 - 64 7% 2 q ) + 55s®) 


+ 


4 

9 


In s 0 (—3 + 5 ?t i 2 q — 18?t i q 


9?77® + s 0 ( — 5 + 45777g — 9777^ — 3?77®) 


9 ■Sq( 1 T 2777^) T 2777gSQ 2 Sq) , 


from which it is easy to extract the moments ((so — m q ) n ). Similarly, weighting with extra 
factors of Eq yields the radiative correction to the moments A4^ n,m \ for n > 1 and any m. 
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The one-loop radiative corrections to the hadronic energy spectrum, and hence to the 
moments , are considerably more difficult to compute. This is because they receive 

contributions from both virtual graphs and bremsstrahlung graphs, only the sum of which 
is infrared finite. The complete calculation of the radiative corrections to the differential 
energy spectrum dT/d-Eo was computed by Czarnecki, Jezabek and Kuhn |L3] . 

We present the leading perturbative corrections to (Eq) and (Eq) in Appendix A. In the 
limit rh q —> 0 they take the simple form 



= 0) = 

[1381 

7 2 I 


-0.768 


900 

30 ^ J 

IT 

71 


= 0) = 

[2257 

4 2 

(y. s 

-0.250 

cx s 

[3600 

45 J 

71 

71 


(3.9) 


C. Corrections to the Moments 

We now combine the results of the previous subsections to present the full expressions 
for the parton-level moments, including the leading perturbative and power corrections. 
The first two moments of the hadronic invariant mass spectrum are given by 


M (1 ’ o) = 


Q? e 


7r 


91 71 


158 


34 


2873 


-1- m„ -m„ H- m„ H- m'l -m„ 

450 18 9 27 9 9 9 18 9 1350 9 


10 


+ 4 ml + —rrC - -m” + — fh° + 


56 ,io . , . 

—m q ) In m q 


(3.10) 


Ai 
2 m 2 
^2 
2 m 2 


40 . 4 § . 6 20 

— m n - m. H- 

9 9 3 9 3 

13 27 97 

-m 2 — 16 rht + 32m® m® H-mi° — 60m: lnm„ — 12m® lnm 

10 2 9 9 9 2 9 10 9 9 9 9 

3 9 45 7 

- — -?h 2 4“ Shi)) ~ 24m® + —m® ~ -ffig 0 + 12 m q In rh q — 36m® In m q 


M m = — 

7 r 


5 

324 

Ai 


137 

- m n 

450 9 


101 


29 


37 


6341 


-m. H- m„ - m° H - m n + 

36 9 81 9 36 9 18 9 8100 


m 


12 


(3.11) 


2 mf, . 

The first mixed moment is 


10 152 , 14. 8 , 0 56,^, / 

—m„ H- rn„ 4-m„ + 4m„ 4-m: lnm 0 

3 9 27 9 3 9 9 135 9 ) \ 

16 16 . 2 4 0 16 . in 16 . 12 

- 45 + - 16m 0 + Wm q ~ m q + 45™* 


128 


m® In rh n 


Ad (1,1) = — 


7T L100 


Ai 


9 209 

+ 


-m„ 

180 9 


149 . 4 4. 6 49 . 8 1457. 10 

-m„ -|— m„ -m° 4-m„ 

108 9 3 9 36 9 2700 9 


23 . 12 
—m„ 
60 9 


7 . 2 10 . 4 4 ,. 2 23 . 10 4 . 12 . . 

5 m q + ~Q m q + y, + 3 m q + —m q + — m q ) lnm fl 


2 ml 


2 m 2 


23 1. 2 10 . 4 

- m„ + 12 m„ 

90 2 9 9 


. c 13 . 8 27 . in 4. 12 

16m® + —m\ - — m q ° + -m q 

20. 6 


(3.12) 


4-12m* In m q + — rh q In m q 
o 

13 3 3 49 4 

-1—m 2 — 4m„-m® 4-m 70 -m 72 + 12m) In m a — 20m® Inm 0 

30 2 9 9 2 9 10 9 3 9 99 99 





































The first two moments of the hadron energy spectrum are given by 


= 


7 5 5 7 

- -m 2 + 8 m q — 8 m® H— m q - rh q ° + 6In rh q + 6m® In m q 

20 4 4 20 

+ AArriq) — 


7r 


Ai 
2 ml 


1 — 8m 2 n + 8m® — m® — 24m:) In m Q 


2ml 


7mi — 20 mi + 20m® — 7m:) 0 + 24m^ In m 0 — 48m^ In m c 


M(°- 2 ) = 


) i i 

--m 2 — 2 mi + 2?h® H —m . 1 . 0 

L15 5 9 9 q 5 q 


— m 12 — 8?h® In m q 
15 q q 


+ A 2 (m q ) 


Of G 


7T 


Ai 

2m'l 


2ml 


43 3 . 2 . 4 9.o 17. 10 2 . 12 

- m„ + 14m„ — 16m® H— m„ - mz° H— m„ 

90 2 9 q q 2 q 10 9 9 9 

28 

+ 12m^ In rhq H-m® In fh q 

3 

13 21 „ . , 21.o 49. in 10. 19 

- mi + 50m„ — 52?n® H-m„ H- ml 0 - m]r 

30 2 9 9 9 2 9 10 9 3 9 


+ 12mf In rh g + 28m® In m c 


(3.13) 


(3.14) 


where the functions A m (rh q ) are presented in Appendix A. 

Finally, we obtain the leading corrections to the total decay rate by taking the n — m — 0 
moment. Of course, this result is not new; we present it for completeness and because we 
will need it to normalize the moments. We find 


r {b - x q ev) = r 0 


fo(m q ) + -^—ofiirhq, Ai, A 2 ) + A 0 (m q ) 


2ml 


a. 


7T 


(3.15) 


where 


fo(m q ) = 1 — 8 m 2 + 8m® — — 24m^ In fri q , 

/i (rh q , Ai, A 2 ) = Ai (l — 8m 2 + 8m® — m® — 24m(j lnm g ) (3.16) 

+ A 2 (—9 + 24 m 2 — 72m q + 72 m® — 15 m s q — 72?h)j lnm g ) . 


The power correction fi(m q ,Xi,X 2 ) was first obtained in Refs. and the perturbative 

correction Ao(m g ), which we present in Appendix A, was first found in Ref. jnj. It takes a 
simple form when rh q —s► 0, for which 


A 0 (m q = 0) 



(3.17) 
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IV. APPLICATION TO B MESON DECAYS 


The relations p.3| ) allow moments of the physical parameters E B and Sh to be expressed 
in terms of the parton-level moments. For the first two moments of Sh we find 

(s H ) = m 2 + A 2 + ( m 2 B - 2Am B + A 2 + + 3A 2 ) (so - in 2 ) 

+ (2A m B — 2A^ — Ai — 3A 2 ) (E $), 

( s 2 h) = m 4 q + 2A 2 m 2 + 2m|(m 2 + A 2 )(s 0 - m 2 q ) 

+ 2m 2 {2Am B - 2A 2 - Ai - 3A 2 )(E 0 ) 

+ (m B - 4Am| + 6A 2 m| + 2X 1 m B + 6A 2 m|)((s 0 - in 2 ) 2 ) (4.1) 

+ 4A 2 m|(Eo) + 4Am|(£ 0 (so - m 2 )), 

where all expressions are valid to relative order Aq CD /7?t| and to all orders in m q /mb. It 
is straightforward to extend the analysis to higher moments (s B )- Similarly, the leading 
moments of E B are given by 

, „ . t Ai + 3A 2 / t Ai T 3A 2 \ . - . 

<B ' ,> = A_ ^^ + ( mB_A+ ^^) <Bo> - 
(£|> = A 2 + (2Am B - 2A 2 - Ai - 3A 2 )(A’o) (4.2) 

+ (m^ — 2A?n B + A 2 + Ai + 3X 2 ) (Eq) . 

On the right hand side of these expressions appear the parton level moments {E™Sq), which 
are obtained from the quantities by multiplying by the scale factor r 0 /r. 


A. Decay to an up quark 

For m q = 0, such as in the quark decay b —> u£u, we find the simple expressions 


, 91 a s 

(So) = —— + 


13Ai 


3A 2 



100 7 r 


180m 5 


60 m B 


(4.3) 


accurate up to corrections of order a s /m 2 B and 1 /m%. These then yield the physical moments 


(s H ) = m 


B 


91 a. 


+ 


7A 


450 7 r 10m 


B 


227 a s 
630 7i 


+ 


10m B 


(A 2 + A, 


A, 


10 
























(4.4) 


5 a a 604A a s 8 A 2 _ \A 
324 7T + 2025ms tt + 15 m% \ 3 J 


(Eh) = lm B 
( E h) = J7 m B 


137 a s 13A 

1 H-- H- 

630 7r 7ms 

257 a s 13A 

1 H-- H- 

480 7T 4ms 


+ 137 g a \ 12A 2 

V 1170 it ) ' 7m|J ’ 



17 a<A 
780 ^T/ 


+ 


13 

4m| 




B. Decay to a charm quark 


For b —> c decays, we make use of the fact that the charm quark is also heavy to write 
m c /m b as a power series in 1 /ms, l/m». Let us define the spin-averaged meson masses, 


m D = 

rriB = 


ms + 3ms* 
4 

ms + 3ms* 
4 


m c + A-——— 1975 MeV 

2 ms 

m b + A-——h ... ~ 5313 MeV , 

2m b 


(4.5) 


which gives 

m c _ % A / fh D \ A 2 / ms\ Ai / _ m^ 

m b rfiB vhb V ms/ m# ' ms/ 2msms \ rhs 

A A 2 A 

= 0.372 - 0.628 - 0 . 628 ^- + 1 . 16 —, 

ms m^ mjj 

accurate up to corrections of order 1/msm^. This substitution introduces additional 
0(l/ms, 1 /msms) corrections to the parton level moments. [j We find 



A + 0.51+ + 1.14 A2 


f y 

(s 0 - mj) = 0.051— + 0.16 


7T 


((s 0 - ml) 2 ) = 0.0053— + 0.017 

y 7T 


7T ms 

A 


m 


B 


- 0.14 


Ai 


(A 0 ) = 0.489 
(El) = 0.242 


. a. 


7r ms 

A 


m 


2 ’ 


1 + 0.043— - 0.78-(1 - 0.12— 

7r ms 

a s A 

1 + 0.099— - 1.50- 

7T m B 


1 - 0.12 


oi s 

IT 

a 


TT 


(E 0 (s 0 - ml)) = 0.030— + 0.077—— + 0.18^f + 0.53^- , 


a o A 


Ai 


m 


2 ’ 
R 


(4.7) 


0.44+ + 1.96+ + 2.53+ 

mjj m^ 


m| 

A 2 


0.19^2- + 3.64 


Ai 


m 


B 


m 


+ 4.69 


Ao 


R 


m 


B 


it ms 


m 


B 


A 2 


and for the total rate, 


2 In the rest of this section, we will treat ms/ms as 0(1). Thus, by 0(l/ms) we denote corrections 
both of order 1 /ms and 1 /ms- 
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r 


= 0.369 


cd, A 

1 - 1.54—+ 3.35- 

7r ms 


The physical moments are then 


1 - U 


0.051— + 0.23— (1 + 0.43—) + 0.26^ (A 2 + 3.9Ai - 1.2A : 
7r ms V 7r / v 


) — 

7T 


A 2 


Ai 


+ 5.81^ - 5.69—+-7.47 


rriR 




A2 


(4.8) 


m 


B 


CD, A CDs 1 /to 

0.0053— + 0.067-- + 0.065—5- (A 2 - 2.1A 

7T ms vr m- v 


B 


(4.9) 


(s/f - m 2 D ) = m 2 B 

<(4 - ? 4 ) 2 > = 

(H#) = 0.489ms 
{E%) = 0.242m| 


where the corrections to these expressions are of order a s /m 2 B and l/m|. Note that in these 
expansions, there is no hidden dependence on the quark masses; here the coefficients are 
functions only of physical quantities. 

We can also expand our results about the small velocity (SV) limit |2 (J, Aqcd m& — 
m c -C m c < mb. In this limit, only the D and D* states are produced. Expanding in powers 
of 1 — m c and AQco/(cnb — m c ), we find 



CD, 
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( CD A 
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4.3A X 
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1 + 0 . 19 — 

+ 0.33— 7T 

(A 2 
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2.9A 2 ) 
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1 + 
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0.099— 

± 

0.55— 

f i + 0.28—^ 

+ 0.75—^2" 

fA 2 

+ 

3.5Ai 

+ 

2.2A 2 ) 


7 r 


m B 

V 7 T J 

m B 

V 




J 


(s H - m 2 D ) = 77i 2 B (1 - rh c y 


4 a s 

— — + 


4A 


+ 


A 2 — 2Ai 


21 7T nib ~ m c {mb — m c ) 


+ ... 


+ 0{l-m c ) A . (4.10) 


Note that in the SV limit, as expected, the average invariant mass of the final hadron ism^, 
and therefore the D and D* are produced in the ratio 1:3. Furthermore, corrections to the 
average invariant mass due to production of excited states are suppressed by (1 — m c ) 3 . 

Finally, all of these results may be applied to the inclusive decays of the A&, with the 
obvious replacements A —> A\, Ai —> Aia, A 2 —s► 0, where 


m Ah = m b + A a 


A 


1A 


2m/, 


(4.11) 


We also note that, in order to avoid introducing factors of A and Ai from the meson sector into 
the expansion, in Eq. (|4.5|) the spin-averaged meson masses should be replaced by baryon 
masses. Since the uncertainty in m Ab is ±50 MeV |2l|, this introduces large uncertainties 
into the moments of A/, spectra, when written in terms of physical masses. 


V. A LOWER BOUND ON A 


Although the invariant mass spectrum for B —> X c ev has not been measured, we may 
use the recent OPAL measurement |22|] of the branching ratio to the narrow P wave charmed 


mesons, the Hi(2420) and ±>2(2460), to place a lower limit on A. In Ref. [22], the branching 
ratio to these states was estimated to be 34 ± 7%. From Ref. P5| , we take the ratio of D 
to D* production in B —> X c iv , for which several experimental measurements have been 
combined consistently: 
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r (B -*■ D*eu e ) 

Y(B —> Deu e ) + r(R —■> D*eD e ) 


0.65 ±0.06. 


(5.1) 


We estimate the minimum value for the first moment of the invariant mass spectrum by 
taking the la limits of these experimental results. Hence, we take a 27% branching fraction 
to the P wave states, and assume that the rest of the branching fraction is saturated by the 
D and D* in the ratio 0.41:0.59. The minimum value for the first moment (sh — mjj) is 
then 


{s H ~ fh 2 D ) min. ~ 0.27 (2.450 GeV) 2 - (1.975 GeV) 2 ± 0.43 (2.010 GeV) 2 - (1.975 GeV) 2 


0.30 


(1.869 GeV) 2 - (1.975 GeV) s 


(5.2) 


= 0.51 GeV 2 


For the second moment, we will be conservative and neglect the small (and positive) contri¬ 
bution of the ground state doublet. We find 


{(s H - rh 2 D ) 2 ) min . ~ 0.27 x [(2.450 GeV) 2 - (1.975 GeV) 2 
Solving Eq. (|4.9|) for the first moment, we find 


1 2 


= 1.2 GeV 4 . 


A > 


Q/ 

0.41-1.41—-0.07 

7T 


Ai 


0.1 GeV 2 


GeV. 


(5.3) 


(5.4) 


The nonperturbative parameters A and Ai are well-defined only at a given order in pertur¬ 
bation theory |^4|. Our limits apply to these quantities defined at one loop. We will use 
the coupling constant a s (rrib ) = 0.2 in what follows. Since Ai is closely related to minus 
the kinetic energy of the b quark in the B meson, it is expected to be negative. Under this 
assumption, we obtain the lower bound 


A > 340 MeV. 


(5.5) 


This limit corresponds to an upper bound on the b quark pole mass of m^ ole < 4.97 GeV 


In Ref. [^5||, the stringent inequality Ai < 
we would find the more restrictive bound 


-3A 2 


-0.35 GeV was proposed; in such a case 


A > 570 GeV, (5.6) 

corresponding to the upper limit ml° lc < 4.71 GeV. 

If we also use the bound (|5.3|) on the second moment, we may relax the assumptions on 
Ai and obtain correlated limits on A and Ai. These are plotted in Fig. Q. By this method, 
we obtain the lower bound 


A > 410 MeV, 


(5.7) 


independent of Ai. Where the bound on A is saturated, Ai = —0.11 GeV 2 . Our result implies 
the upper limit m^ ole < 4.89 GeV, without any assumption on Ai being made. 

|6 


This approach complements the recent proposal 


that A and Ai be extracted from 
moments of the photon energy spectrum in the rare process B —> AGy/. 
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FIG. 4. Correlated one-loop limits on A and Ai. The shaded region is ruled out by our analysis 
of the first two moments of (sjj — m 2 D ). 


VI. HIGHER LOOPS 


In order to apply our results consistently, it is important to know the scale at which to 
evaluate a s (/i) in the radiative corrections. It has been shown recently [27] that the naive 
choice /i = significantly underestimates the size of the two-loop effects. In particular, 
the prescription of Brodsky, Lepage and Mackenzie (BLM) [14|] suggests that the relevant 
scale for the radiative corrections in B —> X u ev decay is p ~ 0.07 mb , when expressed in 
terms of the b quark pole mass, indicating that two-loop effects are substantial. 

It has also been stressed, however, that the BLM prescription may give a misleadingly 
low scale when relating unphysical quantities In particular, A is related to the 

pole mass of the heavy quark, which is not an observable, and in fact suffers from an 
inherent ambiguity in its definition p3. In this section, we show that although the BLM 


prescription indicates that radiative corrections to the first two moments of the invariant 
mass spectrum for semilcptonic b —> u decay are uncontrolled when expressed in terms of the 
HQET parameter A, they are well behaved when expressed in terms of physical quantities. 

The portion of the two loop correction to Eq. i [l. l| j which is proportional to the QCD 
evolution parameter /?o m ay be determined from the one loop correction, calculated with a 
massive gluon in the final state, using the techniques of Ref. |10] 
computation are given in Appendix B; we End, for rh q = 0, 


Some of the details of the 


m-. 


-(sh) = 


91 a s {rrib) ( 53 2 276043A / a s (mb) 


450 7T 


+ 


V180 


-7T 


^ 0.20^ + 3.15 ' 


7T 



10 ms 


— 0.013 + 0.013 + -7-—— + • • • , 
10 ms 


+ 


7 A 

10 


+ ... 


(6.1) 
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where (3q — 11 — 2n//3 and in the last line we have taken a s (m b ) — 0.2. Clearly the 
perturbation expansion is poorly controlled. In the BLM scale-setting prescription, the 
scale /Ublm of the coupling is chosen such that the two-loop contribution proportional to f3 0 
is absorbed into the one-loop correction. The poor convergence of the series is reflected in 
the low BLM scale for this process: 


A* blm = rn b exp 



53 , 

- 7 r 

180 


276043\ 7 91 ' 
108000/ '450. 


~ 0.03m b ~ 140 MeV. 


( 6 . 2 ) 


However, our expression for (s#) is given in terms of the unphysical parameter A. While 
this is perfectly acceptable as an intermediate step, since we are ultimately interested only 
in relations between observable quantities, it has the effect of making the perturbative 
expansion appear ill-behaved. Instead, let us define the “decay mass” of the b quark, w T h , 
via the charmless semileptonic partial width of the B meson, 

T(B - X u ei> e ) = (6.3) 

The decay mass m\ is a physical observable and is therefore well-defined. It is related to 
the pole mass via the expansion 


m h = m 


pole 


/ 5 _ _ 2 _ 2 \ ®s{m b ) 
V6 15^ ) 7T 


(0.596/3 0 + c) 


'a s (m b )' 


+ ... + 


7T 


Ai - 9A 2 

10m^ 


+ 


(6.4) 


The two-loop term proportional to /3 0 , which one expects to dominate the two loop result, 
was calculated in Ref. [p7 |. The constant c has not been computed. Since m^ olc is not 
well-defined due to renormalon effects, the perturbation series in Eq. ( |6.4|) has a renormalon 
ambiguity at 0(l/m b ). 

Dehning a physical version of the parameter A.[] 

A l = m B ~ ml , 

we have 

2 


A 


A l 


+ -- 


mg mg \6 15 

and Eq. (El becomes 
1 


-7T 


- (0.596ft + c) ' 


7T 


7T 


+ 0(a 3 s ,l/m 2 b ), 


m 


2 (s H ) ^ (0.202 - 0.337 ) as ^ + ( 3.151 - 3 . 752 ) ( a ^ mb ^ 


+ 


B 


7 r 


7 r 


IK 

10 itib 


^ -0.135^4 _ 0 . 6 OI lALH 


+ 


7T 


7T 


IK 

10 ms 


7 A r 

~ -0.0086- 0.0024+ b 


10 7Ub 


(6.5) 

( 6 . 6 ) 


(6.7) 


3 Note that unlike A, viq — is not universal for heavy quarks, and differs in the b and c systems. 
Since it explicitly violates heavy quark symmetry, it is not useful to reformulate HQET in terms 
of this more physical quantity. 
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The perturbation expansion clearly has improved dramatically. The corresponding BLM 
scale is now 


+blm = exp [—(2/9)0.601/0.135] ~ 0.37 rn b , 


( 6 . 8 ) 


which is significantly greater than before. 

It is interesting to note that the cancellation we observe in Eq. (| 6 .7| ) persists at higher 
orders in the bubble sum. Using the techniques of Ref. |3l| we can calculate the n loop bubble 
graph, from which we may extract the coefficient of /3 q a " +1 in the perturbative expansion 
for (sq). Although there is no reason to believe that this is the dominant contribution at 
this order, since there is no /3 0 — 1 ► oo limit of QCD in which the quark and gluon bubble 
graphs dominate, it does give one class of contributions to the n loop graphs which displays 
a factorial divergence at large orders in perturbation theory. 

Using the techniques of Ref. |JT| , the perturbation series in Eq. (p.l|) continues as 


m 


2_ ( S H> = 0.202^ + 3.161 (AfiA 


B 


n 


n 


+940.52 


'a s (m b y 


+ 19347.5 


7T 


+ 51.91 


'a s (m b )' 


7T 


'a s (m b y 


7T 


(6.9) 


7 A 

+ 77T-1 - 

10 m B 


and using the results of Ref. [^9|] for the higher order relation between m b and ml, Eq. (El 
continues as 


-^-(sh) = (0.202 - 0.337) + ( 3.151 - 3 . 752 ) 
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7T 




7T 




+ (51.91 - 50.37) 


' a s (m b y 


7T 


+ (19347.5 - 14424.2) 
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7T 


+148.1 


'a s {m b y 


+ 4923. 
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' a s (m b y 


7 r 


7 A^ 

+ —— + ••• 


7T 


7T 


10 rn B 


~ -0.0086 - 0.0024 + 0.0004 + 0.0026 + 0.0051 + 


IK 

10 m B 


+ 


( 6 . 10 ) 


Note that even at higher orders there is significant cancellation between the two series. This 
is similar to the behaviour observed in a different context in Ref. [^]. The remaining bad 
behaviour presumably reflects the presence of unphysical parameters (such as A 2 and Ai) 
at higher orders in the operator product expansion. Assuming the series is asymptotic, the 
size of the smallest term in the expansion gives a measure of the uncertainty in the sum of 
the series. 

We do not find a similar cancellation for the second moment of s B ■ For rh q = 0 and to 
order 1 /m b , we find 
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( 6 . 11 ) 


( s h) - («o) + 4—— (( E 0 S 0 ) - (sg)) 


m. 


+ ... . 


rri B 


Since (E 0 s 0 ) and (4 are both order a s , there is no /3 0 a 2 term introduced by expressing ( s 2 H ) 
in terms of . However, the naive counting of powers of Ao does not work here, because 
(sq) -C (E 0 so), (so)• Instead, the 0(/3oa 2 ) correction to (sg) is the same order as the 0(a 2 ) 
term introduced by expressing ( s 2 H ) in terms of AJ(. Using the /3 0 a 2 term as an estimate of 
the full two loop correction to (sg) alone, we find, using the same technique as before, 



5 a s (m b ) /277 2 

-i -7T 

324 7T V648 


25511\ 

6075 



ot s (m b y 


+ ... 


71 


( 6 . 12 ) 


and so 
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—(4) = 0.015 

m% 


a s (m b ) 


71 


+ 0.0196/% 


'a s {m b ) s 


( A r 

+ 0.298 — b - -0.48 


a s (m b )\ a s (m b ) 


7T 
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7r 


71 
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(6.13) 


= 0.015+!+ + 0,156 (+5+' 


AF aJm b ) 

— 0.214 —-—+ 


71 


71 


ms vr 


~ 9.5 x 10~ 4 + 6.3 x 10“ 4 - 0.214 Afc + ... . 

ms 7r 


The contribution of the /3q a 2 term to this expression is 7.1 x 10 -4 , and the new 0(a 2 ) terms, 
while of the same order as this one, largely cancel against each other. Since the convergence 
of the perturbation series for (4) still appears to be poor, we may also expect the limits 
on A and Ai which we obtained from ((% — fh 2 D ) 2 ) to be more sensitive to higher order 
perturbative corrections than those obtained from (— Ei 2 D ). 

Finally, note that the appearance of A in (4) is suppressed by a factor of a Sl as is its 
associated renormalon at 0{l/m b ). Since renormalon ambiguities must cancel in relations 
between physical quantities, this means that the large (3oa 2 term in (4) does not correspond 
to a 0{l/m b ) renormalon ambiguity in the perturbation series (|6.12|) . 


VII. SUMMARY AND CONCLUSIONS 

We have used the operator product expansion and the heavy quark limit to compute 
the hadronic energy and invariant mass spectra in semileptonic heavy meson decays. Our 
expressions are complete up to order a s in perturbation theory, and up to order a s /m b and 
1 /ml in the heavy quark expansion. The effects of finite final state quark masses have been 
taken into account, so it is possible to apply our results to the important decay b —> civ. 

Our analysis provides a test of the applicability of the OPE to these decays, and of 
the crucial underlying concept of global duality. Only appropriately weighted integrals of 
the theoretical spectra may be compared meaningfully with experiment, and we focus on 
the leading moments. As an initial application, we used the recent measurement of the B 
branching fraction to excited D mesons to put bounds on the nonperturbative parameters 
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A and Ai. We found A > 410 MeV, which led to a constraint on the b quark pole mass, 
m P ole <- 4.89 GeV. More stringent tests will have to await the availability of more precise 
data. The success or failure of our predictions will determine the confidence with which 
one will trust these theoretical techniques in the extraction of CKM matrix elements from 
semileptonic bottom and charm decays. 

We also investigated the behaviour of the perturbation series at higher order in a s , 
to gain insight into the trustworthiness of the lowest order calculation and the choice of 
renormalization scale ft. We found that when written in terms of the unphysical quantity A, 
the perturbation series for (sh) seems to be quite badly behaved, with a BLM scale /trlm 
too low to be meaningful. However, when we define a more physical “decay mass” m£, and 
through it a physical A£, the perturbation series improves dramatically. The cancellations 
which we find persist to higher order in a s , at least when one includes the leading powers 
of fa. 

We have focused on the application to B decays; however, the BLM analysis suggests that 
the perturbative corrections to (s#) are under control for D decays as well. The extension 
of our results to charm is straightforward. 


ACKNOWLEDGMENTS 

ft is a pleasure to acknowledge helpful conversations with Mark Wise and Lincoln Wolfen- 
stein. This work was initiated at the Weak Interactions Workshop at the Institute for 
Theoretical Physics in Santa Barbara, and we thank the organizers for their gracious and 
efficient hospitality. We are equally grateful to the High Energy Theory group at the Uni¬ 
versity of California, San Diego, for their generous support. M.L. and M.S. also thank 
the Institute for Nuclear Theory at the University of Washington. This work was sup¬ 
ported by the United States Department of Energy under Grant Nos. DE-FG03-90ER40546 
and DE-FG02-91ER40682, by the United States National Science Foundation under Grant 
No. PHY-9404057, and by the Natural Sciences and Engineering Research Council of 
Canada. A.F. also acknowledges the United States National Science Foundation for Na¬ 
tional Young Investigator Award No. PHY-9457916, and the United States Department of 
Energy for Outstanding Junior Investigator Award No. DE-FG02-94ER40869. 


APPENDIX A: THE PARTON LEVEL HADRONIC ENERGY SPECTRUM 


In this appendix we discuss the corrections to the parton level hadronic energy spectrum, 
dPo/dFo- Both the perturbative and the power corrections are somewhat unwieldy; we 
present them here for completeness. 

The power correction may be computed by integrating the doubly differential spec¬ 
trum ( p.2| ) over s 0 . The integral will be nonzero only if E 0 < |(1 + m^), because, as 
discussed in Section II, only in this case does the integration region overlap with the one- 
particle pole at s o = riiq. This is a reflection of the fact that the maximum energy a single 
quark can carry away from the decay is |(1 + rhg). In the presence of additional strongly 
interacting particles such as gluons, the total hadronic energy E 0 can exceed |(1 + rh 2 q ). 
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However, the initial motion of the b quark inside the B meson can produce fluctuations 
of the maximum allowed final quark energy above |(1 + m 2 ). These fluctuations appear in 
the differential rate as singular functions S(E 0 — |(1 + ?h 2 )) and S'(E 0 — |(1 + ?h 2 )), which 
are resummed into a smooth function extending beyond the parton model endpoint. For a 
more detailed discussion of this subject see Refs. [||,|T^,|T^] . 

Including the leading power corrections, then, the expression for the hadronic energy 
spectrum is given by 4 


1 dr 
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16\ E 2 - 


rn„ 


3E 0 — AEq — 2 m 2 + 3E 0 m 
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(1 - m y 5 [E 0 - -(1 + rh 2 ) + 0(a s , 1/mi) 


(Al) 


Integrating this expression with respect to Eq, we find the power corrections (|3.13| ) to the 
moments (E 0 ) and {Eq). 

The expression for the perturbative correction to the hadronic energy spectrum is even 


more cumbersome. For the complete spectrum at finite rh q , we refer the reader to Ref. [13 
As an illustration we present here the perturbative corrections at rh q = 0, separately for 
Eq < | and Eq > 

pert. 
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4 We do not agree with the expression presented in Ref. [12]. 
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FIG. 5. The order a s contribution to the differential energy spectrum (l/To)dr/d£'o, for 
m q = 0, in units of a s /-K. In the region Eq > this is the leading nonzero contribution. The 
logarithmic divergence at Eq = ^ is integrable. 
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+ 16^gLi 2 


( 2E 0 -1 \ 

V 2i 0 ) 


This spectrum is shown in Fig. |5|. The logarithmic divergence as E 0 —■> | from above is 
integrable. The region E 0 > \ receives contributions only from brehmsstrahlung graphs. 
Note that the spectrum falls extremely rapidly with increasing Eq. 

The radiative corrections A m {fh q ) to the moments ( E™) may be obtained by integrating 
the full expressions found in Ref. |T3|. We find 
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The correction A 0 (m q ) to the total rate is equivalent to the result presented in Ref. [|19 


APPENDIX B: BUBBLE GRAPHS 

The n-loop bubble graph contribution to moments of So may be calculated from the one 
loop graph evaluated with a finite gluon mass [[!(],[31]]. I 11 this appendix, we briefly outline 
this calculation using the methods of Ref. [[31]]. Only the bremmstrahlung graphs in Fig. |3| 
contribute to the moments of So for n > 1. We consider the expansion 

dr ~ . fa V +1 

Sr = E^o)^(-?) (Bl) 

d5 0 j=0 \ 7T / 
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where (5q = 11 — 2n//3 and the ellipses denote terms which have fewer powers of /3oa s and 
hence are not obtainable from the bubble graphs. Note that these are not suppressed terms 
in any limit of QCD, although they may be numerically small. The n th moment of So then 
has the expansion 


M (n,0) 


OO 

m f ] Po 

3=0 



(B2) 


where 


= J ds 0 sSdj(so) • 


(S3) 


Dehne d 0 (s 0 , A 2 ) and m^iX 2 ) to be the one-loop corrections calculated with a hnite 
gluon mass A, and A = A /m^,. Then 

4 n) (A 2 )= r W o,A 2 ), (B4) 

J A 2 


and we have |3l 


dj(s 0 ) = 


1 d j 


dui 


sm (mi) r s ° 


u =0 


7 TU JO 


(AV) 


dA 2 


do(s 0 , A 2 ), 


(B5) 


where C is a scheme-dependent constant. I11 the V scheme [|14[], (7 = 0, while in the MS 
scheme, C = —5/3. Eqn. (|B5|) may be written in the form 


dj(s 0 ) ^ ' Ck <7c(so), 

k =0 


(B6) 


where Jk is dehned by 

Jk(s 0) = / dA 2 ln fc (A 2 ) —=-do(so, A 2 ). 
Jo dA 2 


(B7) 


Taking moments of both sides of Eq. (|B6|), we have 


m 


(n) 


3 rl rso /s , /s H /v 

= E C W ds'o Sq / dA 2 ln fc (A 2 )- T -d 0 (s 0 ,A 2 ) 

Jo Jo dA 2 

3 


= 55 °k f dA 2 / ds 0 So ln fe (A 2 ) -4-d 0 (s 0 , A 2 ) 
k=0 Jo J a 2 dA 2 

= E c k [ dA 2 ln A ’(A 2 ) -imJ n) (A 2 ), 

, r, Jo dA 2 


(B8) 


where we have used the fact that d 0 (A 2 , A 2 ) = 0 to move the s 0 integral to the right of the 
A 2 derivative. 

It is straightforward to derive analytic expressions for the moments rriQ l \X 2 ) from the 
graphs in Fig. |3j; however the resulting formulas are lengthy and we will not reproduce them 
here. For j = 1 the integrals in Eq. ( fB^ ) may be performed analytically, giving the O(a 2 /3o ) 
correction to while for j > 1 we performed the integrals numerically to obtain the 

contribution from higher loops in the bubble sum. 
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